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The Gauge-Including Projector Augmented Wave (GIPAW) method, within the Density Func- 
tional Theory (DFT) Generalized Gradient Approximation (GGA) framework, is applied to compute 
solid state NMR parameters for 27 Al in the a, 9, and k aluminium oxide phases and their gibbsite 
and boehmite precursors. The results for well-established crystalline phases compare very well with 
available experimental data and provide confidence in the accuracy of the method. For 7-alumina, 
four structural models proposed in the literature are discussed in terms of their ability to reproduce 
the experimental spectra also reported in the literature. Among the considered models, the Fd3m 
structure proposed by Paglia et al. [Phys. Rev. B 71, 224115 (2005)] shows the best agreement. 
We attempt to link the theoretical NMR parameters to the local geometry. Chemical shifts depend 
on coordination number but no further correlation is found with geometrical parameters. Instead 
our calculations reveal that, within a given coordination number, a linear correlation exists between 
chemical shifts and Born effective charges. 



I. INTRODUCTION 

Aluminium oxide (AI2O3), also known as alumina, is 
one of the most important oxides because of its many 
industrial applications 1 . Corundum (a-A^O^), the most 
stable and common crystalline form of alumina, is the 
final product of the calcination of hydroxides or oxyhy- 
droxides of aluminium at temperatures above 1273 K, 
which can vary according to the precursor. The trans- 
formation from aluminium hydrates to the final oxide is 
not direct and a variety of relatively stable intermediate 
phases can be detected for each combination of starting 
compound and thermal treatment. 

The most commonly used precursors in the synthesis 
of the various alumina phases are gibbsite \^-Al(OH) 3 ], 
bayerite [a-Al(OH) 3 ], and boehmite [j-AlO(OH)]. The 
transformation of gibbsite to corundum can proceed 
through a sequence of hexagonal close packed aluminas 
(x and k) or an alternative sequence in which boehmite 
is initially formed and the a phase is achieved via the 
cubic spinel transition phases (7, 5, and 6) 2 . Bayerite 
can also follow this same sequence through boehmite to 
form corundum 3 or an alternative path transforming to 
r\ and 9 phases . 

Among the transition aluminas, the 7 phase is highly 
valued for industrial applications due to its textural prop- 
erties (surface area, pore volume, pore size), which makes 
it an important material in many industrial processes act- 
ing as an adsorbent, a catalyst and/or catalyst support. 
In petroleum and petrochemical industries 7-alumina 
is used as catalyst support for transition-metal sulfides 
Co(Ni)MoS in hydrotreatment catalysts and metallic al- 
loys in reforming catalysts 6,7 . Due to its low crystallinity 
and the consequent difficulty in characterization, the de- 
bate on the structure of 7-alumina remains open and a 
series of theoretical and experimental works concerning 
this subject have been published along the decades 8-17 . 



Since the transition between the distinct intermediate 
phases is a gradual process, the precise temperature at 
which each phase is obtained with a high degree of pu- 
rity can not be determined by X-ray diffraction (XRD) 
experiments only. Solid-State Nuclear Magnetic Reso- 
nance (SS-NMR) is an important technique for material 
characterization. Long-range order is not a prerequisite 
to distinguish different phases and the knowledge of 27 Al 
NMR data can allow the detection of the onset of phase 
changes during alumina calcination 18 , permitting a dis- 
cussion of the transition mechanisms 19 . A detailed inter- 
pretation of the results remains, however, a challenge. 

The advent of theoretical techniques such as the 
Gauge-Including Projector Augmented Wave (GIPAW) 
method 20 ' 21 enables the ab initio calculation of isotropic 
chemical shielding, <r,- ao , quadrupolar coupling constant, 
Cq, and asymmetry parameter, tjq, in solids. First prin- 
ciples simulations of NMR spectra of structural models 
for transition aluminas are now possible 22,23 and can be 
compared with available experimental data. 

In this paper we apply the GIPAW method to compute 
solid state Magic Angle Spinning (MAS) NMR parame- 
ters for a number of well characterized aluminium oxide 
phases and for their boehmite and gibbsite precursors. 
The results are compared with available experimental 
data to validate the method. In order to contribute to a 
better characterization of 7-alumina, a number of struc- 
tural models, proposed in the literature, are examined 
and their simulated spectra compared to experimental 
ones, thus revealing their adequacy. 

We then examine possible correlations of the predicted 
chemical shifts with local atomic geometry or local elec- 
tronic structure, described through Bader analysis 24 and 
Born dynamical effective charges 25 , finding significant 
correlations. 

The rest of the paper is organized as follows: in Sec. II 
we describe the theoretical methodology and the struc- 
tural models used. In Sec. Ill we present our calculated 
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NMR results, compare them with experiments and dis- 
cuss the resulting correlations. Sec. IV contains our con- 
clusions. 



II. DETAILS OF CALCULATIONS 

A. Electronic structure 

All ab initio calculations in this study were per- 
formed using the codes available within the Quantum 
ESPRESSO distribution 26 , which implements the DFT 27 
framework using a plane waves basis set to expand the 
one-electron wavefunctions of Kohn-Sham equations 28 . 
The effect of exchange-correlation (XC) potential was ex- 
plored by comparing the results of different descriptions 
for this term: Perdew-Burke-Ernzerhof (PBE) general- 
ized gradient approximation 29 and its revision (revPBE) 
by Zhang and Yang 30 . Furthermore, we also considered 
a van der Waals-aware density functional (vdW-DF) 31,32 
recently implemented in Quantum ESPRESSO. Interac- 
tion of valence electrons with nuclei and core electrons 
were treated by the Projector Augmented- Wave (PAW) 33 
method. 

The plane wave kinetic energy cut-off and k-points 
sampling were adjusted to yield less than 1 mRy/atom 
convergence in total energy for all models. A kinetic 
energy cut-off of 45 Ry and expansion of augmentation 
charges up to 220 Ry was sufficient to ensure this crite- 
rion. Integration in the Brillouin zone were determined 
by the Monkhorst-Pack 34 procedure 35 . Both atomic po- 
sitions and cell vectors were fully optimized. NMR chem- 
ical shieldings were converged within less than 1 ppm. 

The Born effective charge tensor Z* a/3 is defined 25 by 
the macroscopic polarization induced, in direction j3 and 
under conditions of zero macroscopic electric field, by a 
zone center phonon displacing atomic sublattice k in di- 
rection a. Effective charges were calculated within den- 
sity functional perturbation theory (DFPT) 36,37 and the 
values presented in the following sections correspond to 
their isotropic component obtained by Z* = Tr[Z* Q/3 /3]. 



B. NMR chemical shifts and quadrupolar 
parameters 

First principles GIPAW calculations 20,21 yield the ab- 
solute chemical shielding tensors for each nucleus, <r(r). 
Isotropic chemical shieldings, <ii SO = Tr[<r/3], are com- 
pared to the experimental isotropic chemical shifts by 
using the standard expression: S iso — a re f — <7i S0 . In this 
work we choose corundum as reference such that 27 Al 
shift of the a phase is aligned to the experimental one at 
ppm. 

The resulting ab initio NMR spectra were obtained 
by using the QuadFit program 38 with the theoretically 
calculated chemical shifts and quadrupolar interaction 



parameters, using the experimental magnetic field inten- 
sity, and normalizing each spectral component to reflect 
the relative number of aluminium types. A Lorentzian 
broadening was added to each spectral feature to obtain 
the best comparison with the experimental lineshape. 



C. Structural models 

Among the oxide phases with well characterized struc- 
tures, we studied the final product a and the transitional 
aluminas 9 and k. 

The structural model for a-alumina, used as a refer- 
ence for all simulated NMR spectra in this work, was pub- 
lished by Ishizawa et al. 39 in an XRD study. The struc- 
ture was reported as a corundum-type, with an hexag- 
onal crystal system and i?3c space group. The crystal- 
lographic cell contains six AI2O3 units in which all alu- 
minium sites are coordinated by six oxygens. 

The 9 phase is present along the 7 — > a transition in 
different dehydration paths and its structure has been 
characterized by Zhou and Snyder 9 , with a monoclinic 
crystal system and C2/m space group. In that work the 
structure was refined from Rietvcld analysis resulting in a 
crystallographic cell with four AI2O3 units, in which half 
of the aluminium atoms are octahedrally (Al oct ) and half 
are tetraedrally (Al te t) coordinated. 

The k phase is one of the intermediate products of 
the dehydration path from gibbsite to corundum. 01- 
livier et al. 40 describes the k phase with an orthorhom- 
bic system and Pna2\ space group. The cell contains six 
AI2O3 units and 25% of the Al 3+ sites are tetrahcdral, 
50% octahedral and 25% in a very distorted octahedral. 

We have also studied the NMR spectra of precursor 
phases gibbsite and bochmitc. For gibbsite we started 
from the structure resolved by Saalfeld and Wedde in a 
single-crystal XRD study 41 . The structure is monoclinic 
with P2\/n space group where the Al 3+ cations are oc- 
tahedrally coordinated by 6 OH~ groups forming double 
layers and occupy two thirds of the octahedral holes in 
alternate layers 42 . The interlayer cohesion is granted by 
hydrogen bonds between these OH~ groups. 

In the also layered boehmite structure, described by 
Christcnscn et al. 43 , each Al 3+ cation is octahedrally co- 
ordinated by 2 OH~ groups and 4 intralayer 2 ~ anions 
in an orthorhombic system with Cmcm space group. As 
in the gibbsite structure, the boehmite double layers in- 
teract with each other via hydrogen bonds, which are 
exclusively interlayer. In this phase, hydrogen bonds are 
organized in chains along [001] direction and different 
bond networks are possible depending on the relative ori- 
entation of neighboring chains. In our calculations we 
considered all possible combinations compatible with a 
2x2 supercell. The structure with the lowest energy was 
found to be the one where nearest neighbor chains are 
antiparallel. However energy differences among various 
combinations are found to be less than 1 mRy/cell which 
is consistent with the experimentally observed disorder 
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in this phase . 

As one of the aims of this work is to contribute to 
the characterization of 7-alumina phase, four different 
7-alumina structural models from published theoretical 
works have been studied. 

The first 7-phase model, J-AI2 03(A), considered here 
was proposed by Gutierrez et al. 45 ' 46 This model, also 
called defect spinel or spinel-like structure, consists of a 
cell with 8 AI2O3 units, in which 37.5% of the cations 
are Al tet and 62.5% are Al oct . Among the 24 O atoms, 
12 are four-fold coordinated (04_/ w) and 12 are three- 
fold coordinated (03-/ o ;d) to aluminium atoms 47 . It is 
important to note that in this model only spinel sites are 
occupied by the Al 3+ cations. 

The next 7-phasc model, r y-Al 2 03(B), used in this 
work was published by Digne et al. 7 and proposed by 
Krokidis et al. 14 . The model has 8 AI2O3 units in the 
cell, 25% of all aluminium atoms are Al tet sites, in a sub- 
lattice of 2 ~ anions. The crystal system is monoclinic, 
but very close to an orthorhombic one with the P2\/m 
space group. In this model cations occupy also non-spinel 
sites. 

Two other 7-phase models were published by 
Paglia et al. 48 ' 49 . The unit cells of these models contain a 
large number of atoms, 32 AI2O3 units, and were gener- 
ated from an extensive search on all structural possibili- 
ties of the J-AI2O3 structures using Fd3>m and IAi/amd 
space groups. In the Fd3m model, ^-AlzO^C), among 
the 64 Al 3+ sites 22 are Al tet , 41 are Al oct and 1 is 
Al pen (five-coordinated), while in the IA\/amd model, 7- 
AI2 03(D), there are 21 Al tet and 43 Al oct sites. In these 
two models, due to the breaking of the local symmetry by 
the variations in cation occupancies and related distor- 
tions in octahedral and tetrahedral sites, the symmetry 
is actually Pi 50 . 



III. RESULTS AND DISCUSSION 

A. Structure optimization and approximations on 
exhange-correlation functional 

Geometry optimization was performed for all phases, 
allowing both the atomic positions and cell vectors to re- 
lax keeping the group symmetry fixed. To investigate the 
effect of approximations on exhange-correlation function- 
als, each geometry optimization was repeated with PBE, 
rev-PBE, and the vdW-DF functionals as described in 
Sec. II A. All XC functionals considered overestimate 
the experimental volume, PBE by about 1.8% on aver- 
age, revPBE by 4.2%. The use of van der Waals func- 
tional was found to have negligible effect, inducing a tiny 
further expansion, on all structures except for gibbsite 
(j-Al(OH)s) [see Table I] where a slight contraction was 
observed. The exception in the case of gibbsite can be 
understood considering its layered and open structure. 

Since ah initio NMR calculations are very sensitive to 
structural details, we performed all NMR calculations 



TABLE I. Volumes of the experimental and optimized cells 
with PBE, rev-PBE, and vdW-DF functionals for the five 
phases with well characterized experimental structures con- 
sidered in this work. 



Structure 


Experiment 


Volume (A 3 ) 
PBE rev-PBE 


vdW-DF 


a-Al 2 3 


254.25° 


261.84 266.30 


268.06 


6-AI2O3 


187.92 b 


192.53 195.54 


196.34 


K-AI2O3 


361.3P 


369.19 375.21 


377.75 


j-AlO(OH) 


261.13 d 


262.61 270.24 


271.65 


•y-Al(OH) 3 


427.98 e 


431.57 449.51 


447.30 



Ref. [39] 6 Ref. [9]. c Ref. [40]. d Ref. [43]. e Ref. [41] 



at optimized positions for each XC functional. In spite 
of the discrepancy in calculated equilibrium volume, the 
average difference in &i SO , \Cq \ and t/q between PBE and 
revPBE were calculated as 1.15 ppm, 0.14 MHz and 0.03, 
respectively. These values were found to have negligible 
effect on total spectra and, except for gibbsite, only the 
PBE spectra will be shown. The numerical values for the 
three XC functionals are reported in Table II. 



B. NMR parameters and spectra for well 
characterized structures 

In this section, we present a detailed comparison of 
our ah initio calculated chemical shifts and quadrupo- 
lar interaction parameters for all the well characterized 
structural phases mentioned in Sec. II C with experimen- 
tal data, as well as with very recent theoretical results 
from the literature. All relevant data are shown in Ta- 
ble II. 



1. a and 9 phases 

Corundum, Q.-AI2O3, shows a single well defined NMR 
peak that we take as reference when comparing spectra 
for other structures. 

For the 8 phase, we compare our results with the ones 
reported in Ref. [19], where a mixture of a and phases 
were examined, and the NMR parameters (see Table II) 
for these two coexisting phases were extracted from the 
experimental spectrum by lineshape fitting. 

The overall good agreement between experimental and 
theoretical spectra can be seen in Fig. 1, upper curves. 
The decomposition in individual contributions is also 
given. The main difference between experiment and the- 
ory is in the asymmetry parameter r]Q for the Al te t site. 
For this peak, it should be noted that a recent theoretical 
work 23 has also reported a value similar to ours. 
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TABLE II. Comparison between chemical shifts and quadrupolar coupling parameters calculated in this work and experimental 
and theoretical data from the literature. 









(ppm) 




Cq (MHz) 








Structure 


previous 


this work 


previous 


this work 


previous 


this work 


Al site 


rLxp. 


in. 


PBE/rev-PBE/vdW-DF 


rjxp. 


111. 


PBE/rev-PBE/vdW-DF 


Exp. 


Th. PBE/rcv-PBE/vdW-DF 


CX-AI2O3 






















U.U {I) 


U.U 


0.0 




9 QQ^ 


2 05/1 98/2 10 


0.00 a 


o 6 


n H/1 /n pi/i /n 

U.U^/ U.U'l/ U.UO 


Q-A1203 






















-3.0 (l.0) a 


-5.9 6 


-4.9/-4.8/-5.0 


3.50 (0.30) a 


3.44 6 


2.91/2.81/2.71 


0.00 (0.10) a 0.18 6 


n 9^ /n /n 9Q 
u.zo/ u.oo/ u.zy 


Altet 


00.0 ^l.UJ 


DZ.O 


62.6/60.2/61.1 


0.4U (^U.lUj 


fi Q/1 ^ 
O.o4 


5.63/5.53/5.37 


0.65 (0.02) a 


0.42 & 


U.^O/U.O/ /u.^z 


K-A1203 




















Al„rt(lV 


o.o c 


0.3 b 


2.0/2.3/2.0 


0.50 c 


-9.98 b 


-9.53/-9.45/-9.41 




0.33 & 


n 9Q /n 90 /n 98 

U.zy/ U.zy/ U.Zo 


Al v+ i{2) 




1.4 


7.1/6.7/6.8 


> 1.0 


O.ZU 


4.43/4.33/4.27 




0.94 6 


n 7fi /n 77 /n as 

U. / 0/ U. f ( / U.Oo 


Al tet (3) f 


68.5 C 


fin it> 
OU. f 


62.6/60.2/61.2 


n 7fi c 
U. / 


c cob 
-O.Oo 


-4 83 /-4 67 /-4 49 


0.3 C 


0.33 fa 


0.33/0.31/0.32 


Al oat {A) f 


5.0 C 


4.4 b 


3.8/3.6/3.5 


0.85 c 


4.51 6 


4 QQ/-4 Qfi/c; f)6 




0.77 6 


0.99/1.00/1.00 


j-Al{OH) 3 


Gibbsite 


















Al oct (l) 


-5.6 d < 9 


2.2 e 


2.2/2.3/2. l h 


4.70±0.20 d 


5.10 e 


-5.30/-4.74/-4.43 


1.00±0.05 d 


0.35 e 


0.32/0.35/0.35 


Al oct (2) 


0.0 d ' 9 


0.0 e 


0.0 h 


2.20±0.20 d 


2.80 e 


2.33/2.06/1.81 


0.75±0.05 d 


0.66 e 


0.81/0.71/0.66 


j-AIO(OH) Boehmite 


















Al oc t 


-1.0 d ' 9 




-0.9/-0.5/-0.2 h 


1.8 - 2.8 d 




2.16/2.08/2.27 


0.5 - 1.0 d 




0.45/0.46/0.67 



a Ref. [19]; 6 Ref. [23]; c Ref. [40]; d Ref. [44]; e Ref. [51]. 



^ Aluminium sites are labeled according to Ref. [40] and the sites symmetry, see Table III. 
9 The experimental 8i 30 (Al 0Ct (2)) = 11.5±0.2 ppm from gibbsite is taken as reference. 
The calculated ai so {Al oc t(2)) = 562.72 ppm from gibbsite is taken as reference. 
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FIG. 1. (color on-line) Comparison between our theoretical 
results and the experimental 27 Al MAS NMR spectra from 
Ref. [19] for a mixture of the a and 8 alumina phases (up- 
per curves). Decomposition of the theoretical spectrum in its 
individual components (lower curves). 



TABLE III. Optimized Al — O distances for the K-AI2O3 struc- 
ture. Additionally, for the experimental 40 and the optimized 
structures, we present the average absolute deviation in the 
distances, Ddist., and in the O — Al — O angles, D ang ., ac- 
cording to equations (1) and (2). 



Al oct (l) Al v+I (2) Al tet (3) Al oct (4) 



Distances (A) 


1.962 


1.991 


1.797 


1.879 




1.961 


2.269 


1.772 


1.959 




1.936 


1.848 


1.756 


1.821 




1.820 


1.829 


1.784 


2.026 




1.917 


2.038 




1.838 




1.968 


1.844 




2.215 


D dist . (A) 










Th. 


0.039 


0.130 


0.013 


0.111 


Exp. 


0.058 


0.171 


0.025 


0.105 


Dang. { ) 










Th. 


5.028 


8.392 


3.572 


7.305 


Exp. 


5.150 


8.590 


5.081 


6.888 
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FIG. 2. (color on-line) Simulated 27 Al MAS NMR spectra for 
K-alumina at 11.7 T (upper graph) and 7.0 T (lower graph). 
Decomposition of the spectra in individual components is also 
given. 

2. k, phase 

The crystal structure of K-alumina has been reported 
in a work by Ollivier et al. 40 , in which 27 Al MAS NMR 
and multiple quantum magic angle spinning (3Q MQ 
MAS) 27 Al NMR experiments were also performed and 
interpreted on the basis of the local structure. Three 
distinct resonance signals were indentified in the spec- 
trum obtained with magnetic field intensity of 7.0 T, two 
of them from Al oct sites. At higher magnetic field of 
11.7 T, further resolution revealed an additional peak for 
this site at high shielding, resulting in the interpretation 
of the complex lineshape for Al oct peak region as the over- 
lap of three distinct peaks: two with high, one with low 
quadrupolar coupling. The high coupling component was 
assigned to the strongly distorted Al oct site, Al(2) in the 
published crystallographic description, and was labeled 
as Aly+i to stress the presence of a strongly elongated 
Al -O bond. 

In order to discuss the local distortion of the alu- 
minium sites we calculated the average absolute devia- 
tion in distances 

Ddist . = T fc^i (1) 

' n 

i=l 

where n = 4 or 6 for Al te t and Al oct sites, respectively. 
The angular deviation of the O — Al — O angles from their 
ideal values in octahedra and tetrahedra are calculated 
as 

^,=E^ (2) 

i— 1 

where k = 6 and j3 re f = 109.47° for the Al tet sites or k 
= 12 and /3 re f — 90° for the Al oct sites. 



As it can be seen in Table III, our calculations confirm 
the order of increasing distortion in octahedral sites sug- 
gested in the experimental analysis and in particular site 
Alv+i(2) is found to be the most distorted. 

The simulated total spectra, reported in Fig. 2, ob- 
tained with our calculated Si SO , Cq and uq, are in very 
good qualitative agreement with the Al MAS NMR ex- 
perimental spectrum of Ref. [40] at both magnetic field 
intensities. 

However, from the theoretical decomposition of the 
composite octahedral peak, it can be seen that the high 
shielding feature around -50 ppm, which was experimen- 
tally observed at high magnetic field intensity and was 
assigned to be Al(2) on the basis of the distortion data, 
is actually part of a bimodal peak with high quadrupolar 
coupling and low asymmetry parameter, which belongs 
to Al(l). Morever, the experimental detection of this 
particular feature only at high field intensity can also 
be understood considering the ab initio spectrum for low 
field intensity, reported in the lower part of Fig. 2. As 
can be seen, at low magnetic field intensity, Al(l) peak 
broadens considerably due to its high quadrupolar cou- 
pling constant, making it diffcicult to observe over the 
background in the experiment. This demonstrates that, 
even though experiments can provide accurate structural 
information, ab initio NMR calculations might be essen- 
tial for an unambiguous peak assignment. 



3. gibbsite and boehmite phases 

The structures and NMR properties of gibbsite and 
boehmite are experimentally well characterized. 44 > 51 ~ 53 . 

As can be seen in Fig. 3, the gibbsite NMR spec- 
trum simulated with the theoretically determined pa- 
rameters agrees very satisfactorily with the experimen- 
tally obtained spectrum by Hill et al. ls . This structure 
is the only one for which we observed a noticeable de- 
pendence of the theoretical spectra on the exchange cor- 
relation functional used, showing improved results when 
using vdW-DF 31 functional (see inset in Fig. 3). In Ta- 
ble II, we present a quantitative comparison between our 
calculated NMR parameters and the experimental results 
reported by Damodaran et al. 44 . Our results for Cq and 
rjQ are in good agreement with both of these experiments 
and can accurately reproduce the asymmetry of the peak. 

In agreement with the previous theoretical work by 
Vyalikh et al. 51 , in Fig. 3 we show that this resonance 
profile can be decomposed as the superposition of two 
distinct peaks with different second-order quadrupolar 
structure, belonging to two distinct Al oct sites in gibbsite. 

As described in Ref. [51] the two sites differ in the OH 
groups surrounding them. This can be seen in Fig. 4 
where a single Al(OH)3 layer is drawn showing that 
among the six OH groups surrounding Al oct (l), two par- 
ticipate in interlayer hydrogen-bonds as donors while the 
remaining four are oriented in-plane and participate to 
interlayer hydrogen-bonds as acceptors. For Al oct (2) the 
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opposite occours. 

We further characterize the two aluminium sites re- 
porting, in Table IV, the 17 O chemical shifts, the Born 
effective charges and the hydrogen-bond connectivity of 
the six oxygen types (labeled in Fig. 4) surrounding them. 
This analysis reveals that the two aluminium sites are 
mostly surrounded by the same types of oxygens and the 
distinction is based on just two different oxygen environ- 
ments (O c and Od), each one neighboring only one type 
of aluminium site. 




Chemical shift (ppm) 



FIG. 3. (color on-line) Comparison between the experimental 
27 Al MAS NMR spectra from Ref. [18] for gibbsite and the 
theoretical spectrum obtained with the vdW-DF functional. 
Decomposition of the spectrum in individual components is 
also given. The effect of different exchange and correlation 
functionals on the simulated spectrum is shown in the inset. 




FIG. 4. (color on-line) View from direction (001) of a single 
Al(OH)3 layer in gibbsite showing the arrangement of the dif- 
ferent types of hydroxyl groups around the aluminium cations 
(large gray atoms). Oxygen (small red atoms) are labeled ac- 
cording to their different environments as referred in Table IV. 

The oxyhydroxide polymorph boehmite, j-AIO(OH), 
was detected by XRD in the experiment by Hill et al. ls 
in the temperature range from 200 to 400° C in a mix- 



TABLE IV. Effective charges, chemichal shifts and hydrogen- 
bonding properties (ID = interlayer donor; IA = interlayer 
acceptor) for the different types of oxygen atoms in gibbsite. 
Labelling according to Fig. 4. 



Site 


z* 




H-Bond 


Oa 


-1.46 


238.2 


ID 


o b 


-1.45 


242.4 


ID 


Oc 


-1.43 


246.2 


ID 


o rf 


-1.41 


229.4 


IA 


o e 


-1.38 


238.4 


IA 


0/ 


-1.36 


228.0 


IA 



ture with gibbsite or the transitional x phase. A direct 
comparison between a simulated NMR spectrum for this 
phase with that experiment is therefore not possible. In 
the experiment by Damodaran et al. AA the lack of high 
resolution in the boehmite NMR spectrum is suggested 
to be due to the disorder in Al positions in the sample 
used, leading to small variations in the isotropic chemi- 
cal shift and quadrupolar couplings. Nevertheless, it was 
possible to extract experimental estimates for Cq and tjq , 
which are in good agreement with our calculated values, 
as presented in Table II. 



C. The 7-A/2O3 phase and its structural models 

As previously mentioned, 7-alumina structure is not 
yet well characterized in the literature and a vast dis- 
cussion about this subject persists. Experimental NMR 
spectra of 7-alumina have been reported in a number 
of studies. In the work by O'Dell et al. 19 , the samples 
obtained with calcination temperatures from 600 °C to 
1000 °C are attributed to mixtures of cubic spinel tran- 
sitional phases, as indicated by the presence of a tetra- 
hedral peak around 60 ppm (taking the a phase as a 
reference). The fivefold-coordinated peak that has been 
observed in this phase has been attributed to inherent 
disorder 19 or surface atom contributions 54 . The MAS 
NMR study of Pecharroman et al. 55 reported a small per- 
centage of AIO5 sites as well, along with an occupation 
of 76.3% and 21% for Al oct and Al te t sites, respectively. 
However, no evidence of a fivefold-coordinated Al peak 
is present in some other experiments, such as the one of 
Ref. [56], obtained from a hydrogenated boehmite pre- 
cursor calcinated at 600 °C for several hours and the one 
of Ref. [18], in which the sample was obtained from the 
calcination of gibbsite at 700 °C through the formation 
of boehmite at 300 °C. 

To gain further insight on the structure of this phase, 
we have calculated the ab initio NMR parameters and 
obtained the simulated spectrum for each one of the four 
7-alumina structural models that have been mentioned 
in Section II C . From the comparison of the calculated 
spectra with the experimental ones obtained for calcina- 
tion conditions associated to 7-alumina, an evaluation of 
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Chemical shift (ppm) Chemical shift (ppm) 



FIG. 5. (color on-line) Comparison of the theoretical spectra (dashed red lines) for the four 7-alumina models considered in 
this study and the experimental 27 A1 MAS NMR spectra (black solid lines) from Ref. [19] (left panel) and Ref. [18] (right 
panel). A Lorentzian broadening of 0.5 KHz was added to each component which were properly normalized to reflect relative 
site abundance. 



the adequacy of these structural models will be possible. 
Since NMR is a sensitive probe of the local structural 
and chemical environment, comparisons based on NMR 
parameters will provide complementary constraints with 
respect to X-ray and neutron diffraction methods used, 
for instance, in the detailed structural search of Ref. [49] 

In Fig. 5 we compare the theoretical spectra with the 
experimental data by O' Dell et al. 19 (left panel, 14.1 T), 
and by Hill et al. ls (right panel, 9.4 T). We also compared 
with the 9.4 T spectrum by Pecharroman et al. 55 (not 
shown) which agrees very well with the one by Hill et al.. 

From Fig. 5, it is clear that the simulated MAS NMR 
spectra for the ^-A^O-^iC) structural model proposed 
by Paglia et al. 49 best reproduces all the experimental 
results considered. The agreement is very satisfactory 
and even better for the more precise spectrum at the 
higher magnetic field [left panel of Fig. 5 ] . 

It may appear natural that the two models that dis- 
play lesser agreement with the experimental spectra are 
the two (models A and B) whose unit cells only contain 
8 AI2O3 formula units and a reduced number of non- 
equivalent aluminium environments, while model C con- 
tains a large number of AI2O3 formula units allowing a 
distribution of NMR parameters, which is more suitable 
for the representation of the broad features observed ex- 
perimentally. Notice however that having a large num- 
ber of non-equivalent aluminium environments may be 
considered a necessary, but not a sufficient condition to 
properly reproduce the experimental NMR spectra. In 
fact both Paglia's structural models (C and D) contain 
64 AI2O3 formula units and reproduce neutron diffrac- 
tion data 49 equally well but only model C (the one asso- 
ciated to Fd3m symmetry) satisfactorily reproduces the 
spectral region of the tetrahedrally coordinated Al atoms 
and even for this model the sharpness of the octahedral- 



TABLE V. Average values and corresponding standard devia- 
tions for 27 Al chemical shifts, absolute quadrupolar coupling, 
asymmetry parameter, and site occupation, decomposed ac- 
cording to Al coordination number, for the four 7-alumina 
models considered in this work. 





Siso (ppm) 


\Cq\ (MHz) 


VQ 


% 


y-Al 2 3 (A) 










AlOe 


-0.6±1.8 


4.90±0.47 


0.86 ± 0.12 


62.5 


AIO4 


52.1±3.0 


6.36±2.65 


0.14 ± 0.08 


37.5 


y-Al 2 3 (B) 










AlOe 


-0.7±4.9 


7.98±3.25 


0.41±0.23 


75 


AIO4 


52.0±0.1 


4.51±0.91 


0.47±0.35 


25 


j-Al 2 3 (C) 










AlOe 


-0.8±4.3 


5.93±2.45 


0.52±0.26 


64.0 


A10 5 


15.9 


7.43 


0.34 


1.6 


AIO4 


59.9±4.7 


7.70±3.00 


0.64±0.25 


34.4 


■y-Al 2 3 (D) 










AlOa 


-0.1±4.2 


6.25±2.38 


0.57±0.23 


67.2 


AIO4 


52.7±6.5 


8.27±2.92 


0.58±0.27 


32.8 



Al peak is not completely satisfactory. Table V shows 
that, on average, model C gives a higher average value 
for 6i SO for the Al te t peak while the other models display 
rather similar values that underestimate the experiment. 

Furthermore, in the J-AI2O3 structural model C, a 
truly AIO5 site is evident, and its calculated Si SO of 15.9 
ppm is easily separated from the peaks of octahedrally (- 
9.4 to 9.4 ppm) and tetraedrally (42.6 to 68.0 ppm) coor- 
dinated sites. Although the small occupation of this site 
in model C cannot reproduce the AlO§ peak observed 
by some experiments, the Si SO is in the experimentally 
measured range, implying that an increase in AIO5 sites 
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due to defects and/or surface effects would explain the 
observed feature. 

It should be noted that in the study by Pechar- 
roman et al. 55 , the complex NMR spectrum in the 7- 
alumina region was analyzed as superposition of a small 
number of peaks for which the quadrupolar interaction 
parameters were estimated for tetrahedral (Cq = 4.7 - 4.9 
MHz) and octahedral (Cq = 3.6 - 3.9 MHz) sites. These 
estimated values do not agree with the ones obtained for 
model C. Since model C shows very good agreement for 
the total spectrum, this implies that extracting NMR pa- 
rameters from a complex spectrum, in absence of further 
experimental or theoretical charaterization, is likely an 
unreliable procedure as it can result in widely different 
distributions of NMR parameters for the same spectrum. 



D. Correlation analysis 

In the literature there have been several attempts to 
correlate features in NMR spectra with local structural 
details such as coordination numbers 22 , local atomic 
bondlcngth and angular distortion 57 , shortest bond- 
length 58 , Mulliken charge population 22 etc. 

Our calculations confirm the dependence of chemical 
shifts on aluminium coordination numbers, such that 4- 
, 5- and 6-coordinated aluminium sites show chemical 
shifts in well separated ranges: the lower the coordination 
number the larger is the chemical shift. 

However, the internal distribution within a given co- 
ordination number docs not correlate with local geomet- 
ric descriptors such as local bondlength, angular distor- 
tions 57 , or shortest bondlcngth 58 . To explain this we 
performed a simple test in the non-spinel model B for 
7- alumina by displacing an oxygen atom while keeping 
the rest fixed and we observed that the NMR parame- 
ters were affected in a region extending up to the third 
coordination shell, thus demonstrating the sensitivity of 
NMR to non-local structural details. 

The only properties that we found to display a sig- 
nificant correlation with the NMR shieldings were alu- 
minium Bader charges and Born dynamical effective 
charges. For both charges the correlation is linear and 
structure independent as demostrated by the linear fits 
obeyed equally well by all phases, as shown in Fig. 6. 

It should be mentioned that Bader charges and Born 
effective charges correlate very strongly with each other 
so that no additional information can be gained by con- 
sidering the two charges together. This is because 7- 
alumina is a crystal with a high degree of ionicity, where 
these two quantities are strongly related. 

Furthemore, we can understand the correlation of 



NMR shieldings with Born effective charges because they 
both originate from local electronic susceptibility, one be- 
ing determined from the current induced by the magnetic 
field, the other measuring the charge flow associated to 
vibrational motion. 

As an hystorical curiosity, we mention that propos- 
ing a correlation between NMR chemical shifts and effec- 
tive charges is not a complete novelty since in a couple 
of papers 59 ' 60 in the late seventies a correlation was re- 
ported to exist for binary semiconductors among Szigeti's 
effective charges and 27 Al and 31 P NMR chemical shifts. 

IV. CONCLUSIONS 

We investigated from first principles the 27 Al NMR 
properties of several well characterized crystalline phases 
of AI2O3 and of two of its calcination precursor phases, 
obtaining very good agreement with available experimen- 
tal results. New insight for the peak assignments in the 
spectra were proposed for some structures. 

This gave confidence in the theoretical approach and 
allowed us to address the open problem of the struc- 
tural characterization of the technologically important 
7-alumina phase by comparing the experimental spec- 
tra with the theoretical predictions calculated for four 
structural models recently appeared in the literature: our 
study supports the model structure with Fdim symme- 
try proposed by Paglia et al. in Ref. [49] as the one 
that best reproduces the NMR experimental results in 
the bulk. 

Calculations confirm that chemical shifts strongly de- 
pend on coordination number. Moreover, within a given 
coordination number, a linear correlation exists between 
chemical shifts and Born effective charges or Bader 
charges. 
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